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Abstract. Bounded-size rules(BSR) are dynamic random graph processes which incorpo- 
rate limited choice along with randomness in the evolution of the system. Typically one starts 
with the empty graph and at each stage two edges are chosen uniformly at random. One of 
the two edges is then placed into the system according to a decision rule based on the sizes of 
the components containing the four vertices. For bounded-size rules, all components of size 
greater than some fixed K > 1 are accorded the same treatment. Writing BSR(f) for the 
state of the system with [nt/2\ edges, Spencer and Wormald [24] proved that for such rules, 
there exists a (rule dependent) critical time t c such that when t < t c the size of the largest 
component is of order logn while for t > t c , the size of the largest component is of order n. 
In this work we obtain upper bounds (that hold with high probability) of order n 27 log 4 n, 
on the size of the largest component, at time instants t„ — t c — n ', where 7 £ (0,1/4). 
This result for the barely subcritical regime forms a key ingredient in the study undertaken 
in [4] , of the asymptotic dynamic behavior of the process describing the vector of component 
sizes and associated complexity of the components for such random graph models in the 
critical scaling window. The proof uses a coupling of BSR processes with a certain family 
of inhomogeneous random graphs with vertices in the type space R+ x X>([0, 00) : No) where 
X>([0, 00) : No) is the Skorohod D-space of functions that are right continuous and have left 
limits, with values in the space of nonnegative integers No, equipped with the usual Skoro- 
hod topology. The coupling construction also gives an alternative characterization (than the 
usual explosion time of the susceptibility function) of the critical time t c for the emergence 
of the giant component in terms of the operator norm of integral operators on certain L 2 
spaces. 



1. Introduction 

The classical Erdos-Renyi random graph can be thought of as a dynamic random graph 
process on the vertex set [n] := {1,2, ... ,n}, where one starts with the empty graph n 
(the graph with n vertices and no edges) and at each discrete time step chooses an edge 
uniformly at random and places it in the configuration. Denote by ER (n) (i) for the state of 
the graph obtained after \nt/2\ steps. For any graph G, denote Cj(G) for the i-th largest 
component, and |Cj(G)| for its size (number of vertices). Classical results ([7,11,12]) say that 
for fixed t < 1, the size of the largest component |Ci(ER (n) (i))| is O(logn) while for t > 1, 
|Ci(ER (n, (t))| ~ f{t)n. Here f{t) > is the survival probability of an associated supercritical 
branching process. For t > 1, the size of the second largest component |C2(ER (n) (t))| = 
O(logn). The largest component is often referred to as as the giant component. 
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There have been several works aimed at understanding the nature and emergence of this 
giant as t transitions from below to above t c = 1 ([2, 13]). In recent years, motivated by 
a question of Achlioptas, there has been a significant interest in investigating more general 
dynamical random graph models. The driving theme has been to understand the role of 
limited choice along with randomness in the evolution of the network, in particular the time 
and nature of emergence of the giant component. The simplest such model that has been 
rigorously analyzed, referred to as the Bohman-Frieze process, can be described as follows. 
Start with the empty graph at time t = 0. At each discrete time step, choose two edges e±, e<i 
uniformly at random amongst all pairs of ordered edges. Place edge e\ = (t>i,t>2) if both end 
points V\,V2 are isolated vertices (components of size one), otherwise use edge ei- 

Despite this conceptually simple modification of the standard Erdos-Renyi random graph 
process, a rigorous understanding of this process turns out to be non-trivial. Write BF (n) (i) 
for the state of the system when we have placed \nt/2\ edges. Bohman and Frieze ([5,6]) 
showed that there exists a time t > 1 and e G (0, 1) such that the size of the largest component 
|Ci(BF (n) (i))| = o{n e ). Thus this simple modification delays the time of emergence of a giant 
component. 

In [24] , these results were substantially refined and extended to the context of all bounded- 
size rules (BSR) which we now describe. 

The bounded-size K-rule process {BSR (n) (i)}i>o. Fix K > 0, this will be a parameter 
in the construction of the process. Bounded-size rules treat all components of size greater than 
K in an identical fashion. Let SIk = {1, 2, . . . , K, zu}. Conceptually w represents components 
of size greater than K. Given a graph G and a vertex v £ G, write C V (G) for the component 
that contains v. Define 



For a quadruple of (not necessarily distinct) vertices ^1,^2,^3,^4, write v for the ordered 
quadruple v = (vt, v 2 , v 3 , v 4 ). Let c G (v) = (cg(«i), cg(v 2 ), cg(v 3 ), cq^))- Fix F C Sl A K . 
The set F will be another parameter in the construction of the process. The -F-bounded-size 
rule(i ? -BSR) is defined as follows: 

(a) At time k = start with the empty graph BSRq := n on [n] vertices. 

(b) For k > 0, having constructed the graph BSR^ , construct BSR^Tj as follows. Choose 
4 vertices v = (vi,V2,V3,V4) uniformly at random amongst all n 4 possible quadruples and 



let c k (v) = c BS n k (v). If c k (v) e F then BSR^ = BSRj™' U { Vl ,v 2 ) else BSR<' 
BSR^U(«3,«4). 



We are interested in the dynamics of the rescaled process BSR (n) (i) = BSRj^/^ . Mathemat- 
ically it is more convenient to work with a formulation in which edges are added at Poissonian 
time instants rather than at fixed discrete times. More precisely, we will consider the follow- 
ing random graph process (denoted once more as BSR (n) (t)). For every quadruple of vertices 
v = (vi, V2, v 3 , V4) £ [n] 4 , let Vtf be a Poisson process with rate l/2n 3 , independent between 
quadruples. Note that this implies that the rate of creation of edges is n 4 x l/2n 3 = n/2. 
Thus we have sped up time by a factor n/2 as in the above discrete time construction. Start 
with BSR (n) (0) = n . For any t > 0, at which there is a point in 7-V for a quadruple v G [n] 4 , 
define 




(1.1) 





otherwise 



(1.2) 



where ct~(v) 



PBSR(t-)(*0- 
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Two examples of such processes are Erdos-Renyi process (here K = 0, Ox = {w} and 
F = {(w,va,w,w)}) and Bohman-Frieze process (here K = 1, £Ik = an d i 7 = 

{(1, 1, j3, j'4) : j'3, J4 G ^if})- Spencer and Wormald[24] showed that every bounded-size rules 
exhibits a phase transition similar to the Erdos-Renyi random graph process. More precisely, 
write C- (t) for the z'-th largest component in BSR (n) (t), and |Q for the size of this 
component. Define the susceptibility function 

00 

S 2 (t) = Y,\Cl n) (t)\ 2 . (1.3) 

i=l 

Then [24] proves the following result. 

Theorem 1.1 (Theorem 1.1, [24]). Fix F G T/ien /or i/ie random graph process associ- 
ated with the F-BSR, there exists deterministic monotonically increasing function S2{t) and 
a critical time t c such that limm c S2(t) = 00 and 

S2 (t) p 

> S2{t) as n — )• 00, for all t G [0, t c ). 

n 

For fixed t < t c , \C{ n) (t)\ = O(logn) while for t > t c , \C{ n) (t)\ = Q P (n). 

Here we use o, O, in the usual manner. Given a sequence of random variables {£ n }n>i 
and a function f(n), we say £ n = 0(f) if there is a constant C such that £ n < C f(n) with 
high probability (whp), and we say £ n = £l(f) if there is a constant C such that £ n > Cf(n) 
whp. Say that £ n = 0(/) if £n = O(f) and £ ra = 0(/). In addition, we say £ ra = o(/) if 
0. 

Thus as i transitions from less than t c to greater than i c , the size of the largest component 
jumps from size O(logn) to a giant component @(n). The aim of this work is to study the 
barely subcritical regime, i.e. to analyze the behavior of the size of the largest component at 
times t = t c — e n where e n — > 0. The following is the main result. 

Theorem 1.2 (Barely subcritical regime). Fix F C and 7 G (0,1/4). Then there 
exists B G (0, 00) such that, 

P ||Cf \t)\ < B^^, Vt < t c - n'A -»• 1, 

as n 00. 

As another consequence of our proofs, we obtain an alternative characterization of the 
critical time for a bounded-size rule given in Theorem 1.3 below. Let X = [0, 00) x T>([0, 00) : 
No) where D([0,oo) : No) is the Skorohod .D-space of functions that are right continuous 
and have left limits with values in the space of nonnegative integers, equipped with the 
usual Skorohod topology. Given a finite measure \i on (X,B(X)) and a measurable map 
k : X x X — > [0,oo) satisfying J XxX k 2 (x, y)/i(dx)^(dy) < 00, define the integral operator 
K : L 2 (X,fi) -> L 2 (X,fi) as 

/C/(x) = / «(x,y)/(y) At (dy),/ G L 2 (A», x G AT. 

We refer to k as a kernel on X x X and /C as the integral operator associated with (k, fi). We 
will show the following result. 

Theorem 1.3 (Characterization of the critical time). Fix F C O^. T/ien i/iere exists 
a collection of kernels {Kt}t>o on X x X and finite measures {^t\t>o on (X,B(X)) such that 
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the integral operators ICt associated with {nf,Pt), t > 0, have the property that the operator 
norms p(t) = \\Kt\\ are continuous and strictly increasing in t. Furthermore, t c is the unique 
time instant such that p(t c ) = 1. 

See Section 4.3 for a precise definition of K t and fj, t . We postpone the discussion of the 
connection between the integral operators in Theorem 1.3 and the BSR processes to Section 
2. Using arguments similar to [17] for the Bohman-Frieze model, one can show that for any 
fixed e > 0, the size of the largest component at time t = t c — e can be lower bounded as 
\C { ?\t)\ > A\ogn/{t c - t) 2 where A is a constant independent of e. Thus up to logarithmic 
factors one expects the upper bound in Theorem 1.2 to be tight. 

Theorem 1.2 plays a central role in the study of the asymptotic dynamic behavior of the 
process describing the vector of component sizes and associated surpluses for BSR processes in 
the critical scaling window and its connections with the augmented multiplicative coalescent 
process. This study is the subject of [4] to which we refer the reader for details. 

1.1. Organization of the paper. The paper is organized as follows. In Section 2 we give a 
discussion of the main result. Section 3 collects some notation used in this work. In Section 
4 we introduce and analyze certain inhomogeneous random graph processes associated with 
the BSR process. Finally, in Section 5 we complete the proofs of Theorems 1.2 and 1.3. 

2. Discussion 

We now give some background, open problems and general discussion of the results in this 
work. 

2.1. Subcritical and supercritical random graphs: There has been a considerable inter- 
est in understanding various properties of random graph models in the barely subcritical and 
supercritical regime, see for example ([15, 16,21]) for various results on complex network mod- 
els such as the configuration model in the subcritical regime, [10, 18] for structural properties 
of such graphs including mixing times of the nearly supercritical Erdos-Renyi random graph, 
[25] for an analysis of the Hamming cube near the critical regime and [8] for an extensive 
analysis of a general class of the inhomogeneous random graphs. In recent years there has 
been a significant effort in understanding a special (non-bounded-size) rule called the Product 
rule ([1]), where one uses the edge that minimizes the product of the components at the end 
points. Simulations in [1] suggest that the nature of the emergence of the giant component 
is different from that observed for rules such as the Erdos-Renyi or Bohman-Frieze process. 
Conceptually such rules tend to be harder to analyze since one needs to keep track not just 
of vertices in components upto size K for K < oo but for all K. There has been recent 
progress in understanding such rules [22]. The subcritical regime for such processes has been 
studied in [23] where it is shown that there exists a critical time t c such that for t > t c , the 
susceptibility function defined as in (1.3) satisfies S2(t)/n — > oo as n — > oo. Furthermore 
there exist functions {/fc( - )}fc>i suc h that for t < t c the proportion of vertices in components 
of size k remains closely concentrated about fk(t) as n — > oo. 

For bounded-size rules, the only known results in the barely subcritical regime are in the 
context of the Bohman-Frieze process in [17] and [3] (also see [14] where scaling exponents 
for susceptibility functions in the Bohman-Frieze process were derived). In [17] it was shown 
that for fixed e, the largest component in the Bohman-Frieze process at time t c — e satisfies 
|Ci(BF (n) (i c — e))| = 0(logn/e 2 ). In [3], upper bounds when e = e n , where e n — > 0, were 
obtained and a result analogous to Theorem 1.2 was shown using the special structure of 
certain differential equations associated with the Bohman-Frieze process. 
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2.2. Optimal scaling and open questions: Note that for a fixed e > 0, Theorem 1.2 
gives an upper bound of -B(log n) 4 /e 2 . One would expect, similar to the special case of the 
Bohman-Frieze process treated in [17], that for fixed e, |Ci(BF (n) (t c - e))| = 9(logn/e 2 ). It 
would be interesting to see if the results in the paper can be refined to prove this result for 
general bounded-size rules. Also note that Theorem 1.2 considers time instants t n = t c — n~ 7 
for 7 G (0, 1/4). A separate analysis that uses this result as a starting point is carried out in 
[4] to treat the component sizes in the critical scaling window i.e. when 7 = 1/3. It would be 
interesting to extend the analysis of the current paper to get refined estimates on the size of 
the largest component when 7 £ [1/4, 1/3). 

2.3. Connection to the discrete time system: We use the continuous time construction 
given in terms of Poisson processes, as opposed to the discrete time construction, for math- 
ematical convenience (see e.g. the various martingale estimates in Section 4). It is easy to 
show that in the continuous time construction, by time t n = t c — n~ 7 , the number of edges 
in the system is of the order nt c — n 1-7 + 0(^/n). Using this and the monotonicity of the 
process it is easy to check that Theorem 1.2 holds for the discrete time version as well. 

3. Notation 

We collect some notation used through the rest of the paper. All unspecified limits are 

p d 

taken as n — > 00. We use — > and — > to denote convergence in probability and in distribution 
respectively. Given a sequence of events {E n } n >i, we say E n occurs with high probability 
(whp) if F{E n } -> 1. 

For a set S and a function g : S — > M fc , we write ||g||oo = Yli=i su PseS bi( s )l> where 
g = (#1, • • • <?fc). For a Polish space S, we denote by BM(5), the space of bounded measurable 
functions on S (equipped with the Borel sigma-field B(S).) For a finite set S, \S\ denotes the 
number of elements in the set. No is the set of nonnegative integers. For ease of notation, we 
shall often suppress the dependence on n and shall write for example BSR(t) = BSR (n) (t). 
Recall the Poisson processes V$ used to construct BSR(-) in Section 1. Let {^Ft\t>G ^ e 
the associated filtration: Ft = a \V${s) : s <t,v E [n] 4 }. We shall often deal with {Ft}- 
semimartingales {J(t)} t> Q of the form 

dJ{t) := a(t)dt + dM(t), (3.1) 

where M is a {Ft} local martingale. We shall denote a = d(J) and M = M(J). For a local 
martingale M(t), we shall write (M,M)(t) for the predictable quadratic variation process 
namely the predictable process of bounded variation such that M(t) 2 - (M, M)(t) is a local 
martingale. 

4. Inhomogeneous random graphs 

Fix K > and a general bounded-size rule F C and recall that {BSR(t)} i>0 denotes 
the continuous time bounded-size rule process started with the empty graph at t = 0. Note 
that K = case corresponds to the Erdos-Renyi random graph process for which results such 
as Theorem 1.2 are well known. Thus, henceforth we shall assume K > 1. We begin in Section 
4.1 by analyzing the proportion of vertices in components of size i for i < K. As shown in [24], 
these converge to a set of deterministic functions which can be characterized as the unique 
solution of a set of differential equations. We will need precise rates of convergence for these 
proportions which we establish in Lemma 4.2. We then study the evolution of components of 
size larger than K in Section 4.2. Finally, we relate the evolution of these components to an 
inhomogeneous random graph (IRG) model in Section 4.3. 



6 



BHAMIDI, BUDHIRAJA, AND WANG 



4.1. Density of vertices in components of size bounded by K. Recall from Section 1, 
(1.1) that c t (v) = c BSR( - t )(v), for v £ [n]. For t > and i £ Qk, define 

Xi(t) = | {v £ [n] : c t (v) = %} | and Xl {t) = X t (t)/n. (4.1) 

Following [24], the first step in analyzing bounded-size rules is understanding the evolution of 
Xj(-) as functions of time as n — > oo. Although [24] proves the convergence of x~i{t) as n — > oo, 
we give here a self contained proof of this convergence with precise rates of convergence that 
will be needed in the proof of Theorem 1.2. 

Note that the BSR process changes values at the occurrence of points in the Poisson pro- 
cesses Vrf, v £ [ra] 4 . We call each such occurrence as a 'round' and call a round redun- 
dant if the added edge in that round joins two vertices in the same component. Note 
that such rounds do not have any effect on component sizes or on the vector x(t) = 
(xi(t) , X2(t) , ■ ■ ■ ,xic(t),x^j(t)). We will in fact observe that such rounds are quite rare. We 
now describe the effect of non-redundant rounds on x(-). For j £ and i £ Qk, write 
A(j;i) for the change AXi(t) := Xi(t) —Xi(t—) at an occurrence time t if the chosen quadru- 
ple v £ [n] 4 satisfies ct-{v) = j and the round is not redundant. It is easy to check (see 
Section 2.1, [24]) that, when j = (j u j 2 , j 3 ,j 4 ) £ F, 

A(j;i) = i ■ (l {jl+j2=i} - l {jl= i } - l{j 2 =i}), for 1 < i < K, 
A(j;w) = l {il+i2=ro} (jil {jl < x} + j2l{j 2 <K}), 
with the convention j\ + J2 = w when the sum of ji, j% is greater than K, and j\ + w = 
w + ji = w for all ji £ VLk- For j = (ji, j2, 33, ji) £ F c one uses the second edge {^3,^4} and 
the expressions for A(j;i) are identical to the above, with (^3,^4) replacing {31,32)- Note that 
the corresponding change in the density Xi(t) = Xj(i)/n is given by A(j;i)/n. For j £ 0,^ 
and t > 0, write 

Q(t;j) :={ve [n] 4 : c t (v) = j} . 

Since each quadruple v £ [n] 4 is selected according to the Poisson process V$ with rate l/2n 3 , 
the above description of the jumps of Xi(-) leads to a semi-martingale decomposition of X{ of 
the form (3.1) with 

d(x i )(t) = J2 E ^^MM^t^WI+E E ^ri{c va (t)^c V4 (t)}, 



j£Fv£Q(t;j) j£F c v€Q(t;j) 

where C v (t) := C„(BSR(i)) denotes the component containing v in BSR(t). 

Define for i £ Qk, the functions F{ : [0, 1]^ +1 — > R mapping the vector 
(xi,x 2 , ...,xk,x w ) £ R K+1 to 



(4.2) 



*?( x ) = 2 E O^ii Via x u + 2 E i ) x h x h x h x 3A- ( 4 - 3 ) 
Note that |A(/;i)| < 2if for all J £ fi^-. Also, 



max 



2 *{CvAt)=C V2 (t)},Yl E l{C, 3 W=C, 4 (t)}^ <n 3 K 

j€Fv€Q(t;j) JeF'#eQ(t;j) 



Thus we have 



2K .„ , 2ivT 2 



|d(xi)(t) - if(x(*))| < ^ • 2Kn 3 = — . (4.4) 

2n 4 n 
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Note that xi(0) = 1 while x;(0) = for other i G tt K . Guided by equations (4.2) - (4.4), [24] 
considered the system of differential equations for x(t) := (xj(t) : j G £Ik) 

a^(t)=if(x(t)), ieQ K , t >0, x(0) = (1,0,..., 0), (4.5) 

and showed the following result. 

Theorem 4.1 (Theorem 2.1, [24]). Equation (4.5) has a unique solution. For all i G £Ik and 

t > 0, Xi(t) > 0. Furthermore X^ieQ K x «(0 = 1 an( ^ limt— >oo ^ro(i) = 1- 

The paper [24] also showed that the functions Xj(t) — )■ x«(t) for each fixed i > 0. We will 
need precise rates of convergence for our proofs for which we establish the following result. 

Lemma 4.2. Fix 5 G (0, 1/2) and T > 0. There exist C\, C2 G (0, 00) such that for all n, 

F ( sup sup \xi(t) - Xi(t)\ >n- 5 ) < d exp (-C 2 n 1 - 25 
\i&n K se[o,T] J v 

Proof. Note that Ff(-) is a Lipchitz function, indeed for x,x G [0, 1]^ +1 , 

\Ff (x) - if (x)| < 4if(iT + l) 4 \ x i ~ £ <l ^ 4 ^(^ + X ) 5 SU P 1^ " £ <l- 

Write L>(i) := sup iG Q K \xi{t) -Xi(t)\ and Mj(i) := M(xi)(t). Using (4.4), we get for all % G 0# 
and t G [0,T], 

/"* IK 2 
\xi(t) - Xi (t)\ < / |iT(x( S )) - if(x( S ))|da + T + |Mi(t)| 

ft 2K 2 

<AK{K + lf D(s)ds + T + |Afi(*)|. 

Jo n 

Taking sup ig Q K on both sides and using GronwalFs lemma we have 

sup D(t)< ( sup sup \M i (t)\ + ^^ ) e ^{K+lfT_ 
*g[o,t] \ien K te[o,r] n / 



Thus, for a suitable di G (0,oo), 



J sup D(t) > n- & \ < VP 



sup |Afi(t)| > din' 6 } . (4.6) 
te[0,T] I 



To complete the proof we will use exponential tail bounds for martingales. From Theorem 5 
in Section 4.13 of [19] we have that, for a square integrable martingale M with M(0) = 0, 
]AM(i)| < c for all t, and (M,M)(T) < Q, a.s., for some c, Q G (0,oo), 

P J sup \M(t)\ > a \ < 2exp i - sup [cxA - Qif>(\)] \ , for all a > 0, 

[o<t<T J I A>0 J 

where ^(A) = eAc ~i~ Ac . Optimizing over A, we get the bound 

D 1 rvr Ml 

(4.7) 



Ol I QIC 

sup |M(t)| > a } < 2 exp <J log 1 + — ) + 

0<t<T 



c V Q 



a Q ( ac 

c-? 1o H 1+ q 



s 
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In our context, note that for any i £ Qki |AAfj(t)| = |Axj(i)| < 2K/n. Also, the total rate 

1 

2K 2 T 



of jumps is bounded by n 4 • J^$. Thus for all i G Qk, the quadratic variation process 



r /2K\ Z n 4 
{ M„M,)(T)<1 (-) XsJ * 



n 

Taking a = din~ d , Q = 2K 2 T/n and c = 2K/n in (4.7) completes the proof. □ 

4.2. Evolution of components of size larger than K. Let BSR*(i) denote the subgraph 
of BSR(i) consisting of components of size greater than K. In this section, we will focus 
on the dynamics and evolution of BSR*(i). Note that BSR*(0) = 0, i.e. a graph with no 
vertices or edges. As time progresses components of size less than K merge and components 
of size greater than K emerge. Three distinct types of events affect the evolution BSR*(i): 

1. Immigration: This occurs when two components of size < K merge into a single com- 
ponent of size > K. We view the resulting component as a new immigrant into BSR*(t). 
Note that the first component to appear in BSR*(t) is an immigrant. 

2. Attachments: This occurs when a component of size < K gets linked to a component 
of size larger than K. The former component enters BSR*(i) via attaching itself to a 
component of size larger than K. 

3. Edge formation: Two distinct components of size larger than K merge into a single 
component via formation of an edge between these components. In this case, the vertex 
set of BSR*(t) remains unchanged. 

We now introduce some functions that describe the rate of occurrence for each of the three 
types of events. For i lt i 2 G Qk, define F£ i2 : [0, 1] K+1 — > M as 



^ii,i 2 ( x ) — 2 I XT x h x h x h x u + x h x h x h x h J ■ (4-8) 

K]eF:{j 1 ,j 2 }={ii,i2} jeF c :{j 3 ,j 4 .}={i 1 ,i2} 



For i\,i2 < K, denote n ■ R^^t) as the rate at which two components of size i\,i2 merge. 
When i\ ^12, this rate is precisely 

(2n 3 )- 4 [ J2 X h(t) X h(t)X j3 (t)X j4 (t) +]T X n (t)X n (t)X 33 (t)X j4 (t)\ := n ■ F? 1>ia (x(t)). 

j£F j£F 

{iiJ2}={ii,i2} 0'3»i4}={ii,ta} 

Thus Ri lt i 2 (t) = Ff j 2 (x(t)). The case i\ = 12 < K is more subtle due to redundant rounds 
linking vertices in the same component. The rate of redundant rounds can be bounded by 
7^3- • Kv? ■ 2 = K, from which it follows that 

\R iA (t)-F*Mt))\<-- 

n 

The case i\ = 12 = vo corresponds to creation of edges in BSR*(t) and n ■ F£, CT (x(t)) is the 
rate of creation of such edges. 

We now give expressions for the rates for the three types of events that govern the evolution 
of BSR*(i). The convention followed for the rest of this section is that for ii,i2 £ ^k- 
ii + i2 = when the sum of is greater than K, and zu + i\ = i\ + w = w for all i\ E O^-. 
I. Immigrating vertices: For 1 < i < K, write 3i(t) := n ■ a*(t) for the rate at which 
components of size K + i immigrate into BSR*(i) at time t. Using the above expressions for 
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the rate of merger of components of various sizes we have 

As before the error is due to redundant rounds which can only occur for i\ = 12 = {K + i) jl 
(and when (K + i)/2 is an integer). Now define functions F® : [0, — > R + , and cij(-) : 
[0,oo) -> [0,oo) by 

i?(x) = J2 ^.faW' = f?(x(t)), (4.10) 

where x(i) is as in (4.5). Then (4.9) says that 

sup |a?(t)-f?(x(t))| <K/n. (4.11) 
te[o,oo) 

Note that for any 5 < 1, the error term in (4.11) is o(n _l5 ). Using this observation along 
with the Lipschitz property of F^ j 2 , we have from Lemma 4.2 that for any fixed T > and 
5 < 1/2, 

P( sup sup \a*(t) -Oi(t)\ > n~ s ) < Ci exp(-CW~ 25 ). (4.12) 

l<i<Kse[0,T] 

The constants C\, C 2 here may be different from those in Lemma 4.2, however for notational 
ease we use the same symbols. 

II. Attachments: Fix 1 < i < K and a vertex v contained in a component in BSR*(i). 
Let, for i < K, c*(t) denote the rate at which a component of size i attaches itself to the 
component of v through an edge connecting the former component to v. This rate can 
be calculated as follows. First note that the total rate of merger between a component of 
size i with a component in BSR*(t) is n ■ F?(x(f)). Since there are X w {t) vertices in 
BSR*(t) each of which has the same probability of being the vertex through which this 
attachment event happens, the rate at which a component of size i attaches to v is given by 
nF? ra (x(t))/A ro (t) = Ff m (5t(t)) / x^(t) . Since is a factor of F^ ro (x), c*(t) is a polynomial 
in x(i). Define the functions F t c : [0, l] K+l -> R and q (•) : R + -)■ M + as 

i^ c (x) = ^(x)/x ro , ci(t) = if (x(t)). (4.13) 

Then c*(t) = F?(x(t)). Once again using Lemma 4.2 we get for any 5 < 1/2 and T > 0, 

P( sup sup |c*(t) — Ci(*)| > n _<5 ) < Ciexp(-C 2 n 1 - 25 ). (4.14) 

i<i<is'se[o,T] 

III. Edge formation: Note that the rate of creation of an edge between vertices in BSR*(t) 
is nF^ ro (x(i)). Since such an edge is equally likely to be between any of the X^(t) pairs of 
vertices in BSR*(i), we have that the rate of creation of an edge between specified vertices 
{vi,v 2 } with vi,v 2 G BSR*(t) is b*(t)/n where b*(t) = F£ >w (x(t))/a4(*). Define F b : 
[0, l] K+1 -> R and &(•) : M+ -)• M+ as 

F fe (x) = i^ )W (x)/a4 and b(t) = F b ( X (t)). (4.15) 

Once more it is clear that F b (~x.) is a polynomial and furthermore b*(t) = F b (5t(t)), so by 
Lemma 4.2, for any 5 < 1/2 and T > 0, 

P( sup \b*(t) - b(t)\ > n- & ) < Ci exp(-C 2 n 1 - 2 ' 5 ). (4.16) 

se[o,T] 



K 
< —. 

n 



(4.9) 
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Write a(i) := {a-i(t)} 1<i<K and c(t) := {ci(t)} 1<i<K . We refer to (a, 6, c) as rate functions. 
In the proposition below we collect some properties of these rate functions. These properties 
are easy consequences of Theorem 4.1. 

Proposition 4.3. (a) For all 1 < % < K and t > 0, b(t),ai(t),a(t) > 0. 

(b) We have 

K K 

Halloo := sup VWi) < 1/2, HcH^ := sup V"cj(i) < 1/2, \\b\loc := sup6(t) < 1/2. 

t>0 ^ t>0 ^ t>0 

(c) lim^ 6(t) = 1/2. 

Proof: Part(a) follows from Theorem 4.1 and the definition of the functions. For (b) 
observe that 

4 

52<H(t) = £i?(x(t)) < ^ E x h(t)*32 (*)*is (t)x U (t) - 2 

Statements on ||c||oo; II^IIcxd follow similarly. 

For (c), note that i 7 ^ (x) > x%/2 since when all the four vertices selected are from 
components of size greater than K, two components of size greater than K will surely be 
linked. From Theorem 4.1 lim^oo x w {b) = 1 and thus limsup^oo b(t) > x^(t)/2. The result 
now follows on combining this with (b). □ 

4.3. Connection to inhomogeneous random graphs. In this section, we describe the 
inhomogeneous random graph (IRG) models that have been studied extensively in [8], and 
then approximate BSR*(t) by a special class of such models. We will in fact use a variation of 
the models in [8] which uses a suitable weight function to measure the volume of a component. 
We begin by defining the basic ingredients in this model. Let X be a Polish space, equipped 
with the Borel cr-field B{X). We shall sometimes refer to this as the type space. Let \i be 
a non-atomic finite measure on X which we shall call the type measure on X. A kernel 
will be a symmetric non-negative product measurable function k : X x X — > R and a weight 
function <p '■ X — > M will be a non-negative measurable function on X. We call such a 
quadruple {X,/j,,k, (f)} a basic structure. 

The inhomogeneous random graph with weight function (IRG): Associated with 
a basic structure {X, /i, k, <j)}, the IRG model RG w (^,|i,K^) is a random graph described 
as follows: 

(a) Vertices: the vertex set V of this random graph is a Poisson point process on the space 
X with intensity measure n/x. 

(b) Edges: an edge is added between vertices x, y G V with probability lA K y , independent 
across different pairs. This defines the random graph. 

(c) Volume: The volume of a component C of RG (n) (Af, /j>, ft, (f>) is defined as 

vol^C) := (4.17) 

zee 

For the rest of this section we take 

X := [0, oo) x W where W := 2?([0, oo) : N ). (4.18) 

We first describe how, for each t > 0, BSR*(t) can be identified with a random graph 
with vertex set in X. Recall the three types of events governing the evolution of BSR*(t), 
described in Section 4.2. Each component in BSR*(i) contains at least one group of K + i 
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vertices, i = 1, • • • ,K which appeared at instant s < t in BSR*(-), as an immigrant. We 
denote the collection of all such groups as Imm(i). For C £ Imm(i), we denote by tq £ (0, t] 
the instant this immigrant appears. Also, to each C £ Imm(i), we associate a path in 
V([0,oo) : No), denoted as wc, such that wc(s) = for all s < tq\ wc(s) = wc(t) for all 
s £ [t, oo); and for s £ [rc,t], wc(s) = \Ac{s)\, where Ac(s) denotes the component that is 
formed by C and all the attachment components that link to C over the time interval [tq , s] . 
Then {{re, wc) ■ C £ Imm(t)} is a point process on X and forms the vertex set of a random 
graph which we denote by T(t). We form edges between any two vertices (tc,wc), (t c ,w' c ) 
in T(t) if the components Acif) and Ac 1 it) are directly linked by some edge in BSR*(t). 

Define, for t > 0, the weight function fa : X —> [0, oo) as 

fa(x) = fa(s, w) = w(t), x = (s, w) £ X. (4.19) 

Note that by construction there is a one to one correspondence between the components in 
BSR*(t) and the components in T(t). For a component Co in BSR*(t), denote by Ic the 
corresponding component in T(t). Then note that 

|C |=vol^(J Co ). (4.20) 

We will now describe inhomogeneous random graph models that approximate T(t) (and 
hence BSR*(i)). Given a set of nonnegative continuous bounded functions a = {aj}i<xi 
j3 and 7 = {7i}i<i<x 011 [0>°°) we construct, for each t > 0, type measures ^(a,/3,7) 
and kernels Kt(a,/3,~f) on X as follows. For i = 1,---K and s £ [0, 00), denote by v s ^ 
the probability law on P([s, 00) : No) of the Markov process {^(^)}rg[s,oo) with infinitesimal 
generator 

K 

(A(u)f)(k) = J2^i( u )(f( k + j)-f( k )), /GBM(No) (4.21) 
i=i 

and initial condition w(s) = K + i. In words, this is a pure jump Markov process which starts 
at time s at state K + i and then at any time instant u > s, has jumps of size j at rate Jj(u). 
Denote by v s ^ the probability law on P([0, 00) : No) of the stochastic process {^(^)} r g[o,oo)i 
defined as 

w(r) = w(r) for r > s , w(r) = otherwise. (4.22) 
Now define the finite measure /jbt(a, f3,~f) = fit as 

f /(x)d^(x) = V f ai {u) ( [ f(u,w)v u ,i(dw)) du, f £ BM(X). (4.23) 
J x ~l Jo \JW / 

Next, define the kernel Kt(ot, /3, 7) = Kt on X x X as 

«t(x,y) = K t ((s,w), (r,w)) = / w(u)w(u)/3(u)du, x = (s,w),y = (r,w) £ X. (4.24) 

J 

With the above choice of /xj, fi^ and with weight function fa as in (4.19) we now construct the 
random graph HG (n) (X, /j, t , K t , fa) which we denote by RG (n) (o:,/3,7)(i). We will refer to 
the set of functions (a,f5,-y) as above, as rate functions. These rate functions will typically 
arise as small perturbations of the functions (a, b, c), thus in view of Proposition 4.3(b) it 
will suffice to consider (a,(3,j) such that maxjUaHoo, ||/3||oo, H7II00} < 1- Throughout this 
work we will assume that all rate functions (and their perturbations) satisfy this bound. 
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The following key result says that for large n, T(t) is suitably close to RG(a, b, c)(i), where 
(a, b, c) are the rate functions introduced below (4.16). In order to state the result precisely, 
we extend the notion of a "subgraph" to the setting with type space X and weight function 
4>. For i = 1,2, consider graphs Gj, with finite vertex set Vj C X and edge set £j. We say Gi 
is a subgraph of G2, and write Gi C G2, if there exists a one to one mapping ^ : V\ — > V2 
such that 

(i) <f)(x) < 0(*(x)), for all x G Vi. 

(ii) {xi,x 2 } G £1 implies {$(xi),$(x 2 )} G «? 2 . 

Lemma 4.4. Fix (5 G (0, 1/2) and Ze< e n = n~ s , n > 1. Define, for t > 0, t/ie set of functions 
a~(t) := {(a;(t) — e n ) V 0} KKif , a + (t) := {a;(i) + £n}i<j<i<: a™d similarly c~(i),c+(f) and 
b~(t),b + (t). Define the inhomogeneous random graphs (IRG) with the above rate functions by 

RG~(i) :=RG(a-,&-,0(i), RG+(t) := RG(a + , b + , c + )(t). 

Then for every T > there exist 0^,0^ G (0,oo), such that for all t G [0, T], there is a 
coupling o/RG~(t), RG + (t) andT{t) such that, 

P{RG~(t) C T(t) C RG+(i)} > 1 - Cgexpj-^n 1 - 25 } . 

Proof: The coupling between the three graphs is done in a manner such that T(t) is 
obtained by a suitable thinning of vertices and edges in RG + (i) and RG~(i) is obtained by 
a thinning of T(t). We will only provide details of the first thinning step. We first construct 
the vertex sets V + and V* in RG + (i) and T(t) respectively. 

Let V + be a Poisson point process on X with intensity n/if, where := /if(a + , 6 + , c + ). 
For a fixed realization of V + , denote by (xf, ■ ■ ■ ,x~^), the points in V + , with xf = (sf,w^~) 
and < sf < • • • s~t < t. Write w + = (wf , • • • w~t). We now construct vertices in the 
corresponding realization of T(t) (denoted as {x%, ■ ■ ■ xjv }), along with the realizations of 
Xi(s), i G {Ik, < s < t, which then defines (a* (s) , b* (s) , c* (s)) , < s < t, j = 1, • • • K , 
as functions of x(s) = (xj(s))j g n x as in Section 4.2. For that, we will construct functions 
Wj : [0,t] -)• N , 1 < j < N and Xi : [0,t] -)• [0,1], i G Q,k- We will only describe the 
construction of Wj,Xi until the first time instant s G (0, t], when the property 

a*(s) < at(s), b*(s) < b + (s), c*{s) < c+(s) for all 1 < k < K (4.25) 

is violated. Denote a for the first time that (4.25) is violated with a taken to be t if the 
property holds for all s G [0,t]. Subsequent to that time instant the construction can be 
done in any fashion that yields the correct probability law for T(t). For simplicity, we assume 
henceforth that a = t. After obtaining the functions Wj,Xi, we set x\ = (r*,t(;*), where r* is 
the first jump instant of Wi (taken to be +00 if there are no jumps) and w* G D([0, 00) : No) 
is defined as w*(s) = Wj(s)lr 0jt i(s) + Wi(t)l^ t oo ^(s). The vertex set V* is then defined as 

V* = {xi, • • • x No } = {x* : t* < t, i = 1, • • • N}. 

We now give the construction of w(s) = (wi(s), ■ ■ ■ wn(s)) and x(s) for s < t. Denote by 
{U}i=i, = to < ti < *2 < ■■■tM < t, the collection of all time instants of jumps of {wf}jL 1 
before time t. Denote by the coordinate of w + that has a jump at time t/., and denote the 
corresponding jump size by jk- We set w(0) = 0, Xj(0) = for i 7^ 1 and xi(0) = 1. The 
construction proceeds recursively over the time intervals (£fc-i,tfc], k = I, - ■ ■ M + 1, where 
tjst+i = t. Suppose that (w(s),x(s)) have been defined for s G [0, for some k > 1. We 
now define these functions over the interval (tk-i,tk]- 
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Step 1: s G (t k -\,t k ). Set w(s) = w(tjt_i). The values of x(s) over the interval will be 
given as a realization of a jump process, for which jumps at time instant s occur at rates 
n- Ri lt i 2 (s), G {1, • ■ ■ K}, ii + i 2 < K, where the function Ri 1: i 2 (s), given as a function of 
x(s) is defined as in Section 4.2. A jump at time instant s, corresponding to the pair (11,12) 
as above, changes the values of x as: 

x h (s) = x h (s-) - — , x i2 (s) = x h (s-) - — , x il+i2 (s) = x h+i2 (s-) + %X n , if i x / i 2 , 
n n n 

— , x il+i2 (s) = ^(s-) H 

n n 

Remaining X{ stay unchanged. The values of a*(s), b*(s), c*(s) are determined accordingly. 



x h (s) = x h (s-) — , x il+i2 (s) = x 2il (s-) + — , if i\ = i 2 . 



Step 2: s = t k . Recall that wf k (t k ) — wf k (t k —) = j k . We define Wi(t k ) = Wi(t k —) for all 
i 7^ ifc. The values of Wi k (t k ) and x(ifc) are determined as follows. 

Case 1: wf k {t k —) = 0. In this case K + 1 < j k < 2i"T and is the first jump of wf k . Define 

for 1 < I < K, Q k (l) '■= X^!=i Ri,jk-i(tk—)i where by definition R^y = if i' > K. Note that 
Q k (K) = a*j k (t k —). We set Q k (0) = 0. The values of Wi k (t k ) and 5t(t k ) are now determined 
according to the realization of an independent Uniform [0, 1] random variable u k as follows. 

• If u k > Q k {K)/a jk (t k -), define (w ik (t k ),x.(t k )) = (w ik (t k -),x(t k -)). 

• Otherwise, suppose 1 < l k < K is such that Q k (l k — 1) < u k < Q k (l k ). Then define 
Wi k {t k ) = j k , Xvj{t k ) = x^(t k -) + 3 -± and 

xi k (t k ) = xi k (t k -) x jk _i k (t k ) = x jk _i k (t k -) - J ' fc - — , ifl k ^j k -l k , 
xi k (t k ) = xi k {t k -) - — , tfh=jk-h- 

The value of all other 2j processes at stay the same as their values at t k —. 
Case 2: wf (t k —) / 0. In this case 1 < j k < K. Once again, the values of Wi k (t k ) and x(t k ) 
are determined according to the realization of an independent Uniform [0, 1] random variable 
u k as follows. 

• If u k > I+5! fc _v+ define ('*)>*(**)) = (™ ifc (t fe -),x(i fc -)). 

• Otherwise, 

w ik (t k ) = w ik (t k -) + j k , x jk (t k ) = x jk (t k -) - — , X m (t k ) = X m (tk~) + —, 

and the value of all other Xi processes stay the same as their value at t k — . 

This completes the construction of (w(s),x(s)) for s £ (t k -i,t k ] and thus by this recursive 
procedure and our earlier discussion we obtain the vertex set 

V* = {x u ■ ■ ■ x No } = {x* : t* < t, i = 1, • • • N}, 

which will be used to construct T(t). 

Having constructed vertex sets V + and V*, we now construct edges. For this we take 
realizations of independent Uniform [0, 1] random variables {«ij}i<i<j<oo and construct edge 

between vertices xf and xf in V + if 

1 j 



1 



rt 

Ui,j < - j b + (s)tv+(s)w+(s)ds. 
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This completes the construction of RG + (i). Finally construct an edge between x* and x*j if 
both vertices are in V* and 

Ui,j < 1 — exp | J b* (s)u>i(s)wj(s)ds 

This completes the construction of T(t). By construction T(t) C RG(a + , b + , c + )(t) on the 
set a = t. Also, from (4.12), (4.16) and (4.14) it follows that F(a < t) < C 3 exp {-C/^} 
for suitable constant 0^,0^. The result follows. □ 
The following is an immediate corollary of Lemma 4.4. 

Corollary 4.5. Fix T > 0. Then with C 3 ,C 4 G (0,oo) and , for t G [0,T], a coupling of 
RG~(t), RG+(t) and T(t) as in Lemma 4.4: 

p{vol^(Xf G -(t)) < vol*(Zf(t)) < vol^ t (Xf G+ (t))} 

> 1 - C 3 exp(-C 4 n 1 ~ 25 ), (4.26) 

where X\ (t) denotes the component in T(t) with the largest volume with respect to the weight 
function fy, and Xf- G (t), If- G+ (t) are defined similarly. 

5. Proof of the main results 

In this section, we will complete the proof of Theorems 1.2 and 1.3. Proof of Theorem 
1.3 is given in Section 5.4 while proof of Theorem 1.2 is given in Section 5.5. Recall that 
Lemma 4.4 says that BSR* can be approximated by RG(a, b, c). Sections 5.2 and 5.3 ana- 
lyze properties of integral operators associated with RG(a, y 5, / y) for a general family of rate 
functions (a,/3,-y). We begin in Section 5.1 by presenting some results about an IRG model 
RG (n) (X , fi, k, (f>) on a general type space X. 

5.1. Preliminary lemmas. In this section, we collect some results about the general IRG 
model ~RG <n) (X, fj,, n, (j>). Let /C be the integral operator associated with (k, //), as defined in 
Section 1. Recall that the operator norm of IC, denoted as ||/C||, is defined as 

II KM 

||/G|| = sup , (5.1) 



where for / G L 2 (A», ||/|| 2 = (J x |/(x)| V(*0) 1/2 . 

Lemma 5.1. Fix (X , fi, K,(f>). Denote the vertex set of ~RG {n) (X , fi, k, <f>) = RG (n) by V n 
which is a rate rifi Poisson point process on X . Let K. be the integral operator associated with 
(k,m). Suppose that ||/C|| < 1 and let A = 1 — \\K,\\. Denote by Xf G the component in RG (n) 
with the largest volume (with respect to the weight function <p). Then the following hold. 

(i) If [|0|| oo < °° an d \\ K \\oo < oo, then for all m G N and D G (0, oo) 

P{«oZ^(Xf G ) >m}< 2nD exp{-CA 2 m} + F(\V n \ > nD), (5.2) 

where C = (81101100(1 + SW^^X)))- 1 . 

(ii) Let for n > 1, A n G B(X) be such that 

g(n) := 8 sup |</>(x)| 1 + 3n(X) sup |k(x, y)| < oo. 
xeA„ V (x,y)eA n xA n J 

Then for all m G N, 

Fivol^l^) >m}< nfi(A c n ) + 2nD exp(-A 2 m/g(n)) + F(\V n \ > nD). 



SUBCRITICAL BOUNDED-SIZE RULES 



15 



Proof: Part (i) has been proved in [3] (see Lemmas 6.12 and 6.13 therein). We now 
prove (ii). Consider the truncated version of RG (n) constructed using the basic structure 
{X,fi,K,(j)} where p, := ^|a„ (i.e. the restriction of \i to A n ), «(x,y) = k(x, y)l An (x)l A „(y) 
and 0(x) = (/>(x)1a„(x). Note that ||k||oo < °o and ||</>||oo < oo- Denote by K, the integral 
operator associated with («,/2). Clearly \\jC\\ < ||/C|| and thus A = 1 — \\}C\\ > A. Consider 
the natural coupling between the truncated and original model by using the vertex set V n := 
V n nA n . Write Xf G for the component with the largest volume in the truncated model. Then 
we have 

P{voV(Jf G ) > m} <F{V n nA^0} + F{V n C A n , vol^(Xf G ) > m} 
=¥{V n nA^0} + P{vol (Xf G ) > m} 

<(1 - exp{-n^(A£ )}) + 2nDexp{-A 2 m/g(n)} +F(\V n \ > nD), 

where the last inequality follows from part (i) and the fact that A < A. □ 
For the proof of the following elementary lemma we refer the reader to Lemma 6.5 in [3]. 

Lemma 5.2. Let k, k' be kernels on a common finite measure space (X, /i), with the associated 
integral operators /C, fC' respectively. Then 

(a) ||/C|| < ||K||2 lAt := [f XxX K 2 (yL,y)fi(dx)ii(dy)) . 

(b) Ifn< k', then\\K\\ <\\K'\\. 

(c) \\JC — JC'\\ < \\IC\\, where IC is the integral operator associated with (\k — 

For the proof of the following lemma we refer the reader to Lemma 6.18 of [3]. 

Lemma 5.3. Let jl,fj, be two finite measures on the space X. Assume jl ^ /U and let g = 
dfi/d^i be the Radon-Nikodym derivative. Let k be a kernel on X x X , and define n as 

«(x,y) := \/g(-x)g(y)R(x,y), xjGl 

Denote by fC [resp. K,] the integral operator on L 2 (X, n) [resp. L 2 (X, jl)] associated with (ft, /i) 
[resp. Then \\IC\\ L 2 ^ = H^H^r^, where 11^11^2^ [resp. 11^11^2^)/ is the norm of the 

operator JC [resp. K] on L 2 (fi) [resp. L 2 (Jl)}. 

We end this section with a lemma drawing a connection between the Yule process and 
the pure jump Markov processes with distribution v s i that arose in the construction of the 
inhomogeneous random graphs RG(q,/3,7), see (4.22). Fix j > 1 and recall that a rate one 
Yule process started at time t = with j individuals is a pure birth Markov process Y{t) 
with Y(0) = j and the rate of going from state i to i + 1 given by := i. Also recall from 
(4.18) that W denotes the Skorohod space W := P([0,oo) : N ). 

Lemma 5.4. Fix 1 < i < K and s > and rate functions a,{3,~f. Let {w(t)}t>o be a pure 
jump Markov process with law v s ^ := v s ^{a, /3,7) as in (4.22). Then 

(i) The process w*(t) := w^/KW^fW^/K can be stochastically dominated by a Yule processes 
Y(-) starting with two particles (i.e. Y(0) =2). 

(ii) Fix t > 0, s £ [0, t] and 1 < i < K. Then we have 

[ [w(t)} 2 iy s4 (dw) < 6K 2 e 2tK ^, 
Jw 

and for any A > we have 

v Sli (w(t) > A) < 2(1 - e -mh\\°o}A/2K _ 
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Proof: Let us first prove (i). Note that under the law u Sj i, the process uu satisfies w(u) = 
for u < s and w(s) = K + i < 2K. Further for times t > s, by (4.21), the jumps of the 
w can be bounded as Aw(t) := w(t) — w(t—) < K at rate at most iu(i)||7||oo- The process 
w*{-) is obtained by rescaling time and space for the process w(-). This is once again a pure 
jump Markov process with jump sizes Aw*(t) < 1 which happen at rate at most one. Further 
w*(0) < 2. This immediately implies that this process is stochastically dominated by a Yule 
process with Y(0) = 2. This completes the proof. 

We now consider (ii). We will use the result in part (i). Note that a Yule process started 
with two individuals at time t = has the same distribution as the sum of two independent 
Yule processes {Yi(i)}t>o and {Y 2 (t)} t >o with Yi(0) = Y 2 (0) = 1. Now fix t > 0, s < t and 
1 < i < K. Let w(-) have distribution v s ^. From (i) we have 

w{t) < d K ■ (Y 1 (tK\\ 1 \\ oc ) + Y 2 {tK\\ 1 \\ oc )). (5.3) 

For simplicity write X\ = Y^tKW'yWoo) and X 2 = Y2(tK\\'~/\\ 00 ). Well known results about 
Yule processes ([20, Chapter 2]) say that the random variables X\ and X 2 have a Geometric 
distribution with p := e~ tK II 7 !! 00 . The first bound in (ii) follows from the Geometric distribu- 
tion and the fact 

I [w(t)] 2 u S)t (dw) < ^ 2 E[(X! + X 2 ) 2 }. 

The second bound in (ii) follows from 

u s ,i({w(t) > A}) < 2F{X 1 > A/2K}. 
This completes the proof. □ 

5.2. Some perturbation estimates for RG(a, b, c). Recall that Lemma 4.4 coupled 
the evolution of T(t) (equivalently BSR*(t)) with two inhomogeneous random graphs 
RG(a + , 6, c + )(i) and RG(a",&,c")(t) which can be considered as perturbations of 
RG(a,6,c)(t). The aim of this section is to understand the effect of such perturbations 
on the associated operator norms. Throughout this section X and <pt are as in (4.18) and 
(4.19), respectively. Given the basic structure {X, m , /i t , c)) t }, t > 0, associated with rate 
functions (a, (3, 7), we denote the norm of the integral operator /Q associated with (m, <pt) as 

pt(a, At)- 

The following proposition which is the main result of this section studies the effect of 
perturbations of (ct,/3, 7) on this norm. For a iT-dimensional vector v = (vi, ■ ■ ■ vk) and 
a scalar 6, v + 9 denotes the vector {v\ + 9, ■ ■ ■ vk + 0) and (v + 9) + denotes the vector 
(( Vi + 9)+,---(vk + 9) + ). 

Proposition 5.5. Fore > let pf = pt(a+e, /3+e, 7+e) and = p t ((a—e) + , (/?— e) + , (7 — 
e) + ), where (a, f3,~f) are rate functions. Assume that max{\\a+e\\ 00 ,\\P+e\\ 00 ,\\'y+e\\ 00 } < 1. 
For every T > 0, there is a C5 S (0, 00) such that for all e > and t 6 [0, T], 

max{\p t - pf\,\Pt- PtD ^ C^yfe- (-loge) 2 . 

Proof of Proposition 5.5 relies on Lemmas 5.6 - 5.10 below, and is given at the end of the 
section. We analyze the effect of perturbation of (3, a and 7 separately in Lemmas 5.6, 5.8) 
and 5.10, respectively. 

Lemma 5.6 (Perturbations of /?). Let (ck,/3,7) be rate functions and /3 £ be be a nonnegative 
function on [0, 00) with supQ< s<00 \(3 £ (s) - f3(s)\ <e. Then 

\pt(ot,P,y) - p t (a,f3 £ ,~f)\ < Ce, 
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where C = 6t 2 K 3 ||o;|| 0O e 2 *ll 7 l!° c . 

Proof: Let (fif,Kt) be the type measure and kernel associated with (a,/?, 7) Note that a 
perturbation in f3 only affects the kernel Kt and not fit. Recall the representation of fit in 
terms of probability measures £ [0, t], i = 1, • • • ,K}. From Lemma 5.4(h) 



I \w{t)] 2 u u ^dw) < 6K 2 e 2tKMac , for all u e [0, t], i = 1, ■ 
Jw 

Denote the kernel obtained by replacing (3 by (3 s in (4.24) by k\. Since 
have from (4.24) that 

ft 



K. 



(5.4) 



< e, we 



F (x, y) | < e / w(u)w(u)du < etw(t)w(t), 
Jo 



fit <8> fit a.e. (x, y) = ((s, w), (r, id)). 

By Lemma 5.2 (a) and (c) we now have 



|Pt(«,/3,7) ~ Pt(a,P e ,l)\ < 



|«t(x,y) - Kf(x,y)| 2 (i^(x)(i^ t (y) 



1/2 



< 



(etw(t)w(t)) dfi t (x)dfi t (y) 

=ety~] / ai(s) / [w(t)] 2 i/ Sji (dw;) 
~l Jo Uw 

<et ■ t\\a\\^ ■ K ■ 6K 2 e 2tKM ™ , 



1/2 



ds 



where the last inequality follows from (5.4). The result follows. □ 
When a or 7 is perturbed, the underlying measure fit changes as well and thus one needs 
to analyze the corresponding Radon-Nikodym derivatives. This is done in the following two 
lemmas. We denote by {Gs}r,< s<00 the canonical filtration on P([0,oo) : No). In the following 
we follow the convention that 0/0 = 1. 

Lemma 5.7. Fix e > and let (a, ,#,7), (a, ^,7) be two sets of rate functions such that for 
all 1 < i < K and s > 0, 

cti(s) - £< di(s) < ai(s), and ji(s) - e < ji(s) < ji(s). 

Fix t > and let fit and ft t be the corresponding type measures on X. For (s,w) £ X 
and j > 1 let r| for the time of the j-th jump of w(-) after time s (fit a.s.). Also write 
A(-u) = w(u) — w(u—) for u > 0. Then there exists eo > such that for all e G (0, So), fi -C fi 
and 

djk f v Saw-a-GO n 7a(t?)(t/) r r K 

-jr( s ' w ) = — x n j:r, s <t- 



x cxp 



w(u) 



i=l 



i=l 



du 



dfi^' 1 ' aA(s)-i<(s) J " 3 7A(r;)lTj 

Proof: For i = 1, • • • K, define finite measures fi l t , f~i\ on X as 

fi\(dudw) = ai(u)v u> i(dw)l[Q tt ]{u)du, fl\{dudw) = dn(u)i> u ^(dw)l[Q^(u)du, 

where v u ^ is defined above (4.23) and v u ^ is defined similarly on replacing ji with 7. We will 
show that 



for all 1 < k < K and s S [0, T], u s f. <C v s ^ and 



dv s ,k 
dv s k 



Li(w) 



(5.5) 
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where 

L s := Uj>i 



,TA(^)(r/) 



l{r;<t} j x exp | - J 



w(u) 



K 



K 



i=l 



du } . (5.6) 



The lemma is an immediate consequence of (5.5) on observing that n\ and \i\ are mutually 

singular when i ^ j, and the relation \i t = Yli=i A*l? tk = Yli=i fit- 

We now show (5.5). Prom the construction of u s f. it follows that, there are counting 
processes {Ni(u)} ue \ s>t -\, i = 1, • • • K, on W such that 



A" 



) = u;(s) + iN%(u), for u G [s, i], a.s. i/ S) fc 



(5.7) 



i=l 



and 
M 



w(r)^i(r)dr under v a ^ is a local martingale for u G [s,t]. 

(5-8) 

From standard results it follows that L l s is a local-martingale and super-martingale (see The- 
orem VI. T2, p. 165 of [9]). In order to show (5.5), it suffices to show that {£"}u£[s,t] is a 
martingale. By a change of variable formula it follows that (see e.g. Theorem A4.T4, p. 337 
of [9]) 

K 



K = i + Y, / L 's 



1=1 



1 dMi(u), v G [s,t\. 



(5.9) 



In order to show L* is a martingale, it then suffices, in view of (5.8), to show that (see e.g. 
Theorem II.T8 in [9]) for all 1 < % < K , 



L u s • |7i(u) - ji(u)\w(u)du 



dv S; k{w) < oo. 



Finally note that L" < e etw ( t > . Using Lemma 5.4(i) and standard estimates for Yule processes, 
it follows that for e sufficiently small sup sg j ^ sup 1<fc<ii - f w w(t)e £tw ^ dv s ^(w) < oo. 

The result follows. □ 

We will now use the above lemma to study the effect of perturbations in cc on pt(ct, /3,'y). 

Lemma 5.8 (Perturbations of a). Fix e > 0. Let (a,(3,~f) be rate functions and let a £ = 
(af , • • • ,a £ K ), where af are continuous nonnegative bounded functions on [0, oo) such that for 
all 1 < i < K and s G [0, T] 

cti(s) — e < af {s) < Oij(s). 

Then for every t > 0, 

\pt(a,f3,j) - pi(o; e ,/3,7)| < Cy/e, 



where C = 4/3^ ■ 6K 2 e 2tK Ml« ■ AtKy/^Z. 

Proof: Let (fj, t ,K t ) be the type measure and kernel associated with (a,/3,7). Also, let /if 
be the type measure associated with (o; £ ,/3,7). By Lemma 5.7, 

g(s,w) := P-(s,w) = °^hg^ for (s,w) G [0,t] x W. 
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Using Lemma 5.2 (c), (a), Lemma 5.3, and the fact that |Kf(x, y)| < t||/3||oow(t)w}(i), pt <8> pH 
a.e. (x,y) = ((s,w),(s,w)), we have 



\pt{a,p,j)- pt(cf,p,i)\ < ( / |v/s(x)s(y) " l| 2 |Kt(x,y)| 2 d/z t (x)d/xt(y) 

IXxX 



1/2 



v i / |y / 5(x)5(y) - l| 2 w 2 (t)u; 2 (t)d / u t (x)^(y) 



1/2 



A' 



< 



tWloodi Yl 



a?(s)t^(«) 



1/2 



i^[o,t]2 1 y a^s^u) 



1 I ai(s)ctj(u)dsdu 



(5.10) 



where 



di = sup sup / \w(t)\ 2 v Sii (dw) < 6K 2 e 2tKlM 

se[0,T] KKxJw 



and the last inequality follows from (5.4). In order to bound (5.10), note that: 



<2v^ (v^iM + ^aj(u)j . 



Plugging the above bound in (5.10) gives the desired result. 



□ 



We will now analyze the effect of perturbations in 7 on pt(a,/3, y). We need the following 
preliminary truncation lemma. 

Lemma 5.9. For every T > 0, there exist Cq,Cj,Aq £ (0,oo) such that for any t G [0,T] 
and rate functions (a,/3,y) the following holds: Let pt, Kt be the type measure and kernel 
associated with (ct,/3,y). Define, for A £ (0, 00), the kernel KA,t as 



«A,i(x,y) = Kt(x,y)l{ w (t)<A,w(t)<A} where x = ( s ' w )iy 



[r,w . 



(5.11) 



Then for all A > Aq, 



p{Kt) - C 6 e- C7A < p(K A , t ) < p(Kf), 



where p(nt) [ resp. p(KA,t)J denotes the norm of the operator associated with (nt,m) [ resp. 
(KA,t,fH)]- 



Proof: The upper bound in the lemma is immediate from Lemma 5.2 (b). We now consider 
the lower bound. For the rest of the proof, we suppress the dependence of K t ,KA,t,pt on t. 
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Note that, from Lemma 5.2 (a,c) 

p(k) - p(k a ) < ( / (k(x, y) - k a (x, y)) 2 a>(x)a>(y) I 
\JxxX / 

^2, 



1/2 

1/2 



<2 / (t\\/3\\ x w(t)w(t)l{w(t) > A}) 2 dfi(x)dfi(y) 

\JXxX 

( K K ^ 

<2t|| ) 5|| 00 {di^y j ai(s)aj(u)dsdu 

<2t||/3|| 0O -t||a|| 0O -7^, (5.12) 
where 

di = / [u;(t)] 2 i/ Si j(du;) / [^)] 2 l{„,( t ) >j4 }iv(cfe;). (5.13) 

By (5.3), w(t) <d K(Xi + X2) where X\,X 2 are independent and identically distributed with 
Geometric p = e _ *^" 7 "°° distribution. 

/ [w{t)fl {w{t)>A} u s ^dw) <K 2 E [(Xi +X 2 ) 2 l {Xl 

+X 2 >A/K}\ 

JW 

=K 2 E [(Xi + X 2 ) 2 (l{x 1 +X 2 >C,X 1 >X 2 } + 1 {X 1 +X 2 >C,X 1 <X 2 })] 
<4fT 2 E [X 2 l{ Xl>A /2K} + X 2 2 1{X 2 >A/2K}] i 

The above quantity can be bounded by 

d2{1 _ e -ZTKM-)A/*K < ^exp -f— ^ 4 

for some constant c^- The result now follows on using the above bound and (5.4) in (5.13) 
and (5.12). □ 

Lemma 5.10 (Perturbations of 7). For every T > 0, f/iere exists Cs £ (0, 00) and £0 € (0, 1) 
such that for all t £ [0, T] and rate functions (a, f3, 7) the following holds: Suppose e £ (0, eo) 
and 7 £ = (7f,-- - 7^-), where 7? are continuous, nonnegative maps on [0, T] saca that for all 
l<i< K 

li(s) -e< j!(s) < 7(5), /or a// s £ [0,T]. 

Taen 

|p t (a,/3,7) -/9 t (a,/3,7 E )| < Cgy^ • (- loge) 2 . 

Proof: Let (//t, Kt) [resp. (yuf , k| )] be the type measure and kernel associated with (at, /3, 7) 
[resp. (q,/?,7 £ )]. By Lemma 5.7, for (s,w) £ [0, t] x W 



-p-( a , w ) = n j > 1 VJ/ i { 







j x exp < 









du 



.i=l i=l 

Denote the right side as L\(w). Then, as in the proof of Lemma 5.7, it follows that 
m(w)} u( z[ Stt ] is a {G u }ue[s,t\ martingale under u S)k for every k = 1, • • • K. Fix A £ (Aq,oo), 
where Aq is as in Lemma 5.9, and let K A t be defined by (5.11). Similarly define K £ At by 
replacing Kt with k\ in (5.11). Denote the operator norm of the integral operators associated 
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with (KA,t,IM) and (k £ a t , pf) by pA,t(ct, (3, 7) and pA,t(ot, (3, 7 e ), respectively. Then, by Lemma 
5.3 and Lemma 5.2 (a,c), 



< 



<tA' 



XxX 



1/2 



(K J 4,t(x,y)) 2 d/i t (x)d^(y) 



K K 

EE 

i=i 3=1 



1/2 



[Q,i]x[0,t] 



WxW 



a>i(s)aj(u) ^Lt(w)Ll(w) - 1 v S)i (dw)v u j{dw) 



(5.14) 



Next, using the martingale property of L*, we have 

V^H^H-i' 2 



V / Z|^)Z*^) - 1 v S)i (dw)v u j(dw) 



2-2 s/L\(w)v s ,i(dw) / y/Li(w)u u ,3(dw) 



w 



w 



<4-2/ V4H^W-2 / 

w)v u Adw) 



w 



(5.15) 



w 



where the inequality on the last line follows on observing that from Jensen's inequality the two 
integrals on the second line are bounded by 1 and using the elementary inequality ai + ai < 
0102 + 1, for a%,a2 £ [0, 1]. We will now estimate the two integrals on the last line of (5.15) 
by using the martingale properties of {£"}«e[s,t] an d the representations (5.7) and (5.9) in 
the proof of Lemma 5.7. For the rest of the proof we write L" as L s (u). By an application of 
Ito's formula, we have that for every k = 1, • • • K , v s \~ a.s. 



^--Ei^KiH- 

s<u<t 

£ / yi^H 



(«) 



'2f(u) 

7i(«) 



^■A 2 ^ 7i («) 



(5.16) 



1 dJVi(u) 



dNi(u) 



7t(«) 



1 d^(u), 



where the second equality follows on observing that for u G (s, i] ^^Lju) 
As in the proof of Lemma 5.7, we can check that for all i, k, 



w 



J y/L s (u) ■ |t| (u) - ji(u)\w(u)du 
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and consequently the last term on the left side of (5.16) is a martingale. Denoting the 
expectation operator corresponding to the probability measure v s ^ on W by E S) fc, we have 

k r , , ,2 

1-E atk [y/L$)] = ? J> S ,A 



i=l 
A' 



1 dN.iu) 



i=i 
1 X 

E 



w(u)'ji(u)du 



8=1 



E 



<- 



<- 



1 

2 



Ke ■ E Stk [^ L s (u)w(u)]du 

'E S!k [L s (u)}E Sik [w\u)}) 1/2 du 



Ke 



<fii. t . (6K V T ™~) 1/2 , 

where the last inequality follows from (5.4). Using the above bound in (5.14) we now have 



\p A>t (a,p,y) - p A>t (a,P,y e )\ <tA< 



t\\ot\ 



2Ket(6K 2 e 2TK ^) 1 / 2 



1/2 



Finally, by Lemma 5.9, we have 

\pt(a, AT) " Pt(c*,f3,~f £ )\ <\pA,tipt,P,i) - PA ,t(<x,Pn £ )\ + 2C 6 e~ C7A 

<dl A 2 e 1 ' 2 + 2C 6 e- c ^, 

where d\ = tA 2 ^^ ■ i||a||oo ■ [2Kt(6K 2 e 2TK ^°°) 1/2 ] 1/2 . The result now follows on taking 
A = — log e in the above display and taking eo sufficiently small (in particular such that 

-log(£ ) >Aq). □ 

Now we combine all the above ingredients to complete the proof of Proposition 5.5. 

Proof of Proposition 5.5: Using Lemma 5.10, 5.6 and 5.8, we get 

\Pt ~ Pt\ <\pt(a + e,(3 + e,y + e) - p t (ot + e,fi + e,7)| 

+ \p t (cx + e,(3 + e,j) - p t {ot + e,P,y)\ + \pt(a + e, /3,y) - p t {a,/3,y)\ 

<C 8 e^ 2 (- loge) 2 + d x e + d 2 e 1 ' 2 , 

where di = 6T 2 K 3 e 2TK and d 2 = 24K 3 T 2 e 2TK . A similar bound holds for \pj - p t \. The 
result follows. □ 

5.3. Effect of time perturbation on pt. Throughout this section we fix rate functions 
(a,(3,y). The aim of this section is to understand the evolution of the operator norm 
Pt(cx,f3,y) as t changes. The main result of the section is Proposition 5.11 which studies 
continuity and differentiability properties of the function p(t) := pt(ot, f3,y), t > 0. 

Proposition 5.11. Suppose that f3(t) > for t > and liminft_>. 00 (3{t) > 0. Then 
(i) p is a continuous strictly increasing function on R + with p(0) = and lim^oo pit) = oo. 
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(ii) There is a unique value t' c = t' c (at, /3,^f) such that p(t' c ) = 1. 

The proof of the proposition relies on the following lemma and is given after the proof of 
the lemma. 

Lemma 5.12. Let < t\ < t% < oo. Then 



Proof: Letting p := pt 2 we have 

\p(h) ~ P(h)\ <( [ («t 2 (x,y) - K tl (x,y)) 2 d^(x)dfi(y)) 

\JXxX J 

< ( I {W\\oow{t 2 )w(t 2 )\t 2 - hlfd^d^y)) 1 

\JXxX J 

<|i 2 - h\ • W\\oo • t 2 \\a\\ 00 • 6K 2 e 2t2K ^, 

where the last inequality once again follows from (5.4). This proves the upper bound. 
Next note that, for p <g> p a.e. (x, y) such that K tl (x,y) ^ 0, we have 

«t a (x,y) _ 1 + w{u)w{u)(3{u)du 
K tA^,y) J* 1 w(u)w(u)(3(u)du 

w^wfa) M tl < u <t 2 • (t 2 - h) 



>1 + 



Thus Kt 2 (x, y) > ( 1 + \t 2 — ti\ ■ mit \ = y 2 ^ti(x, y) which from Lemma 5.2 (b) implies 

p{t 2 ) - pit!) > \t 2 - ti\ — p(ti). 



This completes the proof of the lower bound. □ 
Proof of Proposition 5.11: Since k,q = 0, the property p(0) = is immediate. Also 
Lemma 5.12 shows that p is continuous and strictly increasing. Finally since inf^oo (3{t) > 0, 
there exists 5 > and a t* £ (0,oo) such that for all t > t*, f3(t) > 5. From Lemma 5.12 
we then have, for t > t*, p(t) — p(t*) > ^wmr ■ This proves that p(t) — > oo as t — > oo and 
completes the proof of (i). Part (ii) is immediate from (i). □ 

5.4. Operator norm of RG(a, b, c) and critical time of BSR. In this section we will 
prove Theorem 1.3. Recall that by Lemma 4.4, for any fixed time t, BSR*(£) (more precisely, 
r(i)) can be approximated by perturbations of RG(a, b, c)(t). To estimate the volume of the 
largest component in RG(a, b, c)(t) we will use Lemma 5.1. In order to identify suitable A n 
as in part (ii) of the lemma, we start with the following lemma. 

Lemma 5.13. Let (a,/3,7) be rate functions and let pt be the associated type measure. Fix 
T > 0. Define A G B(X) as A = {(s, w) £ X : w(T) < 1} for I € R+. Then, for every I £ R + 



p t (A c ) < 2TIHU -exp -I 



2K~ 



24 BHAMIDI, BUDHIRAJA, AND WANG 

Proof: Note that 

^(A C ) = V/ ai{u)v Ut i(k c ) < HallooT sup sup z/ Uii (A c ). (5.17) 

u£[0,T] l<i<K 

By (5.3), 

^(A c ) = u u>i ({w : «;(T) > I}) < F(X t + X 2 > l/K) < 2(1 - e -^NI~y/2* 
where Xj are iid with Geom(e _T " 7 "°°) distribution. Using this estimate in (5.17), we have 

«*(A C ) < HallooT • 2(1 - e -TK\hUy/2K _ 
The result follows. □ 



We will now use the above lemma along with Lemma 5.1 to estimate the largest component 
in RG (n) (o:,/3,7)(i). Recall the notation pt(a,f3,~f) from Section 5.2. 

Lemma 5.14. Let (a:,/3,7) be rate functions and denote byT^ G (t) the component with the 
largest volume, with respect to the weight function <j)t, inRG (n) (i) := RG ( "'(a,/3,7)(f). Then, 
for every t > such that pt(a,f3,-f) < 1, there exists A E (0,oo) such that 

P(voZ 0t (T 1 RG (t)) > Alog 4 n) ^ 0, as n oo. 

Proof: We will use Lemma 5.1(h). Define 

A n := {(s,w) £ X : w(i) < Blogn} , 

where B will be chosen appropriately later in the proof. Now consider the function g{n) in 
Lemma 5.1(h) with A n defined as above and (p,4>,K) there replaced by (pt,<ftt, Kt), where 
(pt, K t) is the type measure and kernel associated with (a,f3,~f). Note that 

K t (x,y) = / /3(u)w(u)w(u)du <t\\/3\\ x w(t)w(t) 



and therefore 

g(n) < 8B log n(l + 3p t (X) ■ tW^B 2 log 2 n). (5.18) 
Writing m n = Alog n, the bound in Lemma 5.1(h) then gives 



P(vol 0t (Xf G (t)) > m n ) < np t (A c n ) + 2np t (X)exp{-A z Alog 4 n/g(n) j , (5.19) 

where A = 1 — p t (a,(3,~f) > 0. Using Lemma 5.13 with I = .Blogn gives 

np t (A c n ) < ntHloo • n-Be-™ 00 /™ = (i) ( 5 . 2 0) 
for B > 2Ke T ^ oc . Now fix fi > e T \\i\\°v/ 2K ; and choose A large such that 

n / u i (A')exp^-A 2 Alog 4 n/fif(n)^ ^0 

as n — > 00. The result follows. □ 
Proof of Theorem 1.3: Let, for t > 0, (fj,t,Kt) be the type measure and the kernel 
associated with rate functions (a, 6, c). We will prove Theorem 1.3 with this choice of (ut, Kt). 
From Proposition 5.11 we have that p{t) = pt(a,b,c) is continuous and strictly increasing in 
t and there is a unique t' c £ (0,oo) such that p(t' c ) = 1. It now suffices to show that: (a) For 
t < t' c , \Ci(t)\ (the size of the largest component in BSR*(i)) is 0(log 4 n); and (b) for t > t' c , 
\d(t)\ = O(n). 
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Consider first (a). Fix t < t' c . For 5 > 0, define rate functions (a + ,6 + ,c + ) = (a + S, b + 
5, c + 5). Since p(t) < 1, by Proposition 5.5, we can choose 5 sufficiently small so that 
Pt(a + ,b + ,c + ) < 1. Denote Xf- G+ (t) for the component of the largest volume in RG + (t) := 
RG(a + , b + , c + )(t). From Lemma 5.14 there exists A £ (0, oo) such that 

P(vol 9it (X 1 RG+ (t)) > ,41og 4 n) -»■ 0, as n ^ oo. 

Combining this result with Corollary 4.5 we see that 

P(vol^ t (Xf (t)) > Alog 4 n) ->• 0, as n -> oo, 

where X\{t) is the component with the largest volume in T(t). Part (a) is now immediate from 
the one to one correspondence between the components in T(t) and BSR*(i) (see (4.20)). 

We now consider (b). Fix t > t' c . Then pit) > 1. From Proposition 5.5 we can find 
5 > such that pt(a~, b~, c~) > 1, where (a - , b~, c~) = ((a — 5) + , (b — 5) + , (c — 5) + ). Let 
C^ G (t) be the component in RG~(i) := RG (n) (a~, b~, c~)(t) with the largest number of 
vertices. By Theorem 3.1 of [8], |Cf G ~(t)| = 0(n). Since vol^ (Cf G ~ (t)) > |Cf G "(t)|, we have 
vol^^Xf-* 3 (t)) = U(n), where Xf G (t) is the component with the largest volume in RG~(t). 
Finally, in view of Corollary 4.5, we have the same result with Xf G it) replaced by X\{t) and 
the result follows once more from the one to one correspondence between the components in 
r(i) and BSR*(i). □ 

5.5. Barely subcritical regime for Bounded-size rules. In this section we complete the 
proof of Theorem 1.2. Throughout this section we fix 7 £ (0, 1/4) and let t n = t c — ra -7 . The 
main ingredient in the proof is the following proposition. 

Proposition 5.15. There exist B,C,N £ (0, 00) such that for all n > N and all < t < t n 

C Bdozn) 4 
P{|Cf l) (i)| > fh(n,t)} < -z, where m(n,t) - " ' 



n 2 ' v w (t c -ty ' 

Let us first prove Theorem 1.2 assuming the above proposition. 

Proof of Theorem 1.2: Write r = w£{t > : \C[ n) (t)\ > m(n,t)}, where m(n,t) = 2B ^ t y)* • 
Then 

F{\C[ n) {t)\ > m(n,t) for some t < t n } = P{r < t n }. (5.21) 

Note that 

{r = t}c [J E v,v ' , (5.22) 

v,v' S[n],v^v' 

where, denoting the component in BSR(i) that contains the vertex v £ [n] by Ci n \t) and its 
size by |C^, n) (t)|, 

E?>« ={max{|Cr(t-)|,lC ) (t-)|} < m(n,t);C^ (t-) ^ C ( ;\t-)} 

1^ {|CW(t-)| + \C™(t-)\ > m(n,t)} r|{CW( t ) = C W(t)}. (5.23) 

Note that 

F{p v n \t)\ + \C ( ;\t)\ > m(n,t)} < 2F{\C[ n) (t)\ > m(n,t)/2} (5.24) 

and, on the set, {max {|C£ n (t)\, (t)\\ < m(n,t)}, the rate at which Cy n \t) and C^\t) 
merge can be bounded by 

^•4|CW(t)||C^(t)|n 2 <^^. 
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Combining this observation with (5.22) and (5.24), we have 

P{r < U < Yl + |C«(t)| > m(n,t)} • 2 Jjf^A dt 

<2n 2 /%{|C{ n) (t)| > m(n,t)/2} • 2m ^ n ^ dt 
Jo n 

<4nt c sup {m 2 (n,t)¥{\C[ n) (t)\ > fh(n,t)}} 

t<t n 

=0{n ■ ra 47 (logn) 8 ■ n~ 2 ) = 0(n- 1+47 (log n) 8 ) = o(l), 

where the last line follows from Proposition 5.15 and the fact that 7 < 1/4. Using the above 
estimate in (5.21) we have the result. □ 

We will need the following lemma in the proof of Proposition 5.15. 

Lemma 5.16. Let (a + , b + , c + ) = (a+ 5 n , b + 5 n , c + 5 n ), where 5 n = n -270 and 70 G (7, 1/4). 
Let Pt' + = pt(a + , b + , c + ). TTien i/iere exists Cg, Ao G (0, 00) swc/i i/iat /or a// n > A^o, 

p (n),+ < 1 _ _ Q J Qr a U < t < tn _ 

Proof of Lemma 5.16: From Proposition 5.5, there is a d\ G (0, 00) such that 

P ( t n) ' + < Pt(a, b, c) + din" 70 log 2 n, for all t < t c . 
By Lemma 5.12 and since |0t c (a, b, c) = 1, there exists c?2 G (0, 00) such that 

pt(sL,b,c) < 1 — d%{t c — t), for all t < t n . 
Thus, since 7 < 70, we have for some Nq > 

p W-+ <i- d2 (t c -t) + d in -^(\ognf < 1 - ^(t c - t), 

for all n > A^o and < t < i c — n~ 7 . The result follows. □ 

Proof of Proposition 5.15: Recall the rate functions (a, b, c) introduced in Section 
4.2. Choose 70 G (7,1/4) and let (a + ,o + ,c + ) be as in Lemma 5.16. Fix t < t n and con- 
sider the random graph RG w (a + , b + , c + )(t). From Lemma 4.4, we can couple T(t) and 
RG (n) (a+,6 + ,c+)(f) such that 

P(r(t) C RG ( ' l) (a + ,6 + ,c+)(t)) > 1 - C 3 exp(-C74n 1 - 470 ), for all t G [0,T]. 

Recalling the one to one correspondence between components in BSR*(t) and T(t), and (4.20), 
we have for any m > 1, 

n\C[ n \t)\ >m}< P{vol^(Xf G+ (t)) > m} + C 3 expj-Qn 1 " 470 }, (5.25) 

where Xf ,G+ (t) is the component with the largest volume in RG w (a + , 6 + , c + )(t). From 
Lemma 5.16, there is a Ao > such that A < t n),+ = 1 — /?t(a + , 6 + , c + ) satisfies 

A<™ ),+ > C 9 (t c - t), for all t<t n , n> N . (5.26) 

Using Lemma 5.1 and arguing as in equation (5.19) we have for all t G [0, T] and all m > 1, 

P{vol 0t (Xf G+ (t)) > to} < ndiexp{-(Aj n) ' + ) 2 m/(d 2 log 3 n)} + o^n" 2 , (5.27) 
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where d%, d^,ds are suitable constants. Using (5.26) in (5.27) we get 

P{vol0 t (Xf G+ (t)) >m}< nd 1 exp{-d 4 {t c - t) 2 m/ log 3 n} + d 3 n~ 2 . 

The result now follows on substituting m = m(n,t) = , with B > 3/c?4, in the above 

inequality. □ 
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